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Abstract. 

The realization that GRBs are narrowly beamed implied that the actual rate of GRBs is much larger than the 
observed one. There are 500 unobserved GRBs for each observed one. The lack of a clear trigger makes it hard 
to detect these unobserved GRBs as orphan afterglows. At late time, hundreds or thousands of years after a GRB, 
we expect to observe a GRB remnant (GRBR). These remnants could be distinguished from the more frequent 
SNRs using their different morphology. While SNRs are spherical, GRBRs that arise from a highly collimated 
flow, are expected to be initially nonspehrical. We ask the question for how long can we identify a GRBR among 
the more common SNRs? Using SPH simulations we follow the evolution of a GRBR and calculate the image of 
the remnant produced by bremsstrahlung and by synchrotron emission. We find that the GRBR becomes spherical 
after - 3000yr(£5i/^)i/^ at 7? - 12pc(£5i/^)i/^ where ^51 is the initial energy in units of 10^^ erg and n is the 
surrounding ISM number density in cm' -I We expect 0.5(^51/^)^^^ 

non- spherical GRBs per galaxy. Namely, 
we expect ~ 20 non spherical GRBRs with angular sizes ~ /iarcsec within a distance of lOMpc. These results 
are modified if there is an underlying spherical supernova. In this case the GRBR will remain spherical only for 
~ 150yr(£5i/n)^/^ and the number of non-spherical GRBRs is smaller by a factor of 10 and their size is smaller 
by a factor of 3. 



INTRODUCTION 

A y-ray burst (GRB) that originates within a galactic disk 
deposits ~ 10^^ ergs into the ISM. This results in a blast 
wave whose initial phase produces the afterglow. The late 
phase of the blast wave evolution would result, as noted 
by Chevalier [1] in the context of supernova remnants 
(SNRs), in a cool expanding Hi shell. The shell will 
remain distinct from its surrounding until it has slowed 
down to a velocity of ^ 10km s~^ [2], which should 
happen within 2.3 • lO^yr E'^-^^ where £^51 is the initial 
energy in units of 10^^ erg. 

The observed rate of GRBs is one per ~ 10^ yr per 
galaxy [3]. The implied GRB isotropic energy is of the 
order of 10^^ ergs. These estimates suggested that there 
are a few remnants per galaxy at any given time. As it 
was believed that the GRB explosions were much more 
energetic than SNs, Loeb and Perna [2] suggested that 
GRBRs would form HI supershells. This giant structures 
require much more energy than what a usual SN can 
supply. 

However, the realization that GRBs are beamed [4, 5, 
6, 7, 8] changed both estimate. First the rate of GRBs is 
much higher. Beamed GRBs illuminate only a fraction 
fh of the sky, thus their rate should be higher by a factor 
of f^^. With ft ~ 0.002 [9] we expect several thousand 



GRB remnants per galaxy. On the other hand the energy 
output of each GRB is much smaller [9, 10, 11]. Thus 
they cannot produce the giant HI shells. 

How can we distinguish a GRBR from and SNR and 
for how long? Both GRBs and SNs deposit a compara- 
ble kinetic energy (~ 10^^ erg) into the ISM. The energy 
injection in a GRB is in a form two narrow relativistic 
beams containing ~ 10~^ Mq. A SN deposits this energy 
spherically with ~ lOM©. In both cases the expected late 
evolution is similar. At this late stage both remnants are 
in the Sedov [12] regime where all the kinetic energy 
is in the ejecta and all the mass is in the surrounding 
ISM. A key distinguishing feature unique to GRB rem- 
nants could be their beamed nature. We expect that the 
beamed emission would lead to a distinct double shell 
morphology at intermediate times. The late time behav- 
ior of the GRB remnant is expected to be spherical. To es- 
tablish how many H I shells are GRB remnants we need 
to find out the expected morphology of GRB remnants 
and how long they stay non- spherical and distinguish- 
able from SNRs. Establishing how many of the H I shells 
are GRB remnants would make it possible to directly es- 
timate the local rate of GRBs, determine 8, the efficiency 
of converting the explosion energy into y-rays, and the 
beaming factor [2] . 

We model the intermediate evolution of a beamed 



GRB by two blobs of dense material moving into the 
ISM in opposite directions. We follow numerically the 
hydrodynamic evolution [13]. We find that the morphol- 
ogy depends on the dimensionless ratio between the 
accumulated mass and the initial kinetic energy, ^ = 
McV^o. When ^ - 2.1 x 10^ (at t - 3000yr(£5iA)^/^ 
and R ^ 12pc(£'5i/n)^/^) ,the remnant becomes spheri- 
cal and indistinguishable from a SNR. 

An additional complication arises if the GRB is ac- 
companied by a supernova, as suggested in the CoUap- 
sar model [14, 15, 16]. The supernova produces an un- 
derlying massive spherical Newtonian shell that prop- 
agates outwards. At /i ~ 3000 corresponding to ^ ~ 
150yr(£5i/n)V3 and R 4pc(^5i/^)^/^ this shell will 
catch the non- spherical GRBR and the system will 
quickly become spherical. In this case the number of 
non-spherical GRBRs is smaller by a factor of 10 and 
their size is smaller by a factor of 3. 
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FIGURE 1. A schematic evolution of a GRB during its rela- 
tivitistic phase. From top to bottom: (a) An inner engine accel- 
erates relativistic jets, (b) Collisions within the jets produce the 
observed GRB. (c) External shocks produce the afterglow, (d) 
The jets expand sideways as they slow down. 




THE NUMERICAL SIMULATIONS 



The Model 
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A GRB occurs when a compact 'inner engine' ejects 
two ultra-relativistic beams . Internal collisions within 
these beams leads to the GRB (See a schematic de- 
scription in Fig. 1.). Later external shocks caused by 
collisions with circumstellar matter produce the after- 
glow. The matter slows down during this interaction and 
its bulk Lorenz factor F, decreases. The ejecta stays 
coUimated only until F drops below ~ l/6o? at ap- 
proximately 2.9hr(£5i/m) V^(0o/O.l)^/^ after the GRB 
[18, 4] where 0o is the initial angular width. At this 
time the matter starts expanding sideways causing, for 
an adiabatic evolution, an exponential slowing down 
[18]. The ejecta continues to expand sideways at an 
almost constant radial distance from the source Rq ^ 
0.3 pc El(^n~^^^ until it becomes non-relativistic. At this 
stage, we begin our simulation. 

Without a detailed numerical modeling of the relativis- 
tic phase of the ejecta we have only an approximate de- 
scription of the initial conditions. We expect the angular 
width of the ejecta to be ~ 1 rad and we are constrained 
by the energy conservation: 



Rq ^ 0.3 pc^^^^^^ 



1/3 



(vo/c) 



-2/3 



(1) 



Our initial conditions comprise two identical blobs mov- 
ing at in opposing directions into the ISM. Both 
the blobs and the ISM are modeled by a cold y = 5/3 
ideal gas. The blobs are are denser than the ISM. 

Luckily the intermediate and late evolution of the 
ejecta are insensitive to the initial conditions. Already 
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FIGURE 2. A schematic evolution of a GRBR. From top to 
bottom: (a) Initial conditions around the Newtonian transition, 
(b) Shells collision along the equatorial plane, (c) A late time 
spherical shell. 



in the intermediate stage we are in the Sedov regime, 
the mass is dominated by the "external" ISM gas which 
washes out any variations in the initial conditions of the 
ejecta. Our numerical simulations [13] verified this ex- 
pectation and different initial densities, angular widths 
and shapes of the blobs led to essentialy similar late time 
configurations. 

Our code is based on the Newtonian version of the 
smooth particle hydrodynamics (SPH) code introduced 



in [17]. The code was adapted for the specific problem 
at hand. We have also used the post Newtonian version 
of the code to take account of possible initial relativistic 
effects (with an initial blob velocity of c/3). 

Once we choose the initial velocity. Equation (1) 
leaves us with the freedom of choosing two out of the 
three parameters Eq, Rq and n, the initial energy, distance 
and ISM density respectively. In presenting the results 
we choose £"0 and n. To parameterize the evolution of 
the remnant we utilize the fact that mass scales linearly 
with initial energy and define the dimensionless parame- 
ter jn = Mc^/Eq where M is the accumulated shell mass. 
We define M as all mass with density above 2n. We show 
all subsequent results as functions of jn. The simulations 
begin with ^ ~ 1. Conveniently ju scales linearly with 
time with: t ~ 0.046/i(£5i/^)~^/^yr, as can be seen from 
Fig. 3. 
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FIGURE 4. Equally spaced density contours 
(p = 1.5w, 2«, . . . , 3.5^) Sit ju = 9.5 x 10^, 1.4 x 10^ 2.4 x 
10^ 3.9 X 10^ 6.1 X 10^ 10^, 1.8 x 10^, 3 x 10^, 5 x 10^ (left 
to right, top to bottom) 



FIGURE 3. Time as a function of ju. The linear relation 
between time and jjist ^ 0.046/j{Esi/n)~^^^ yr. 



Results 

As each blob collides with the ISM it produces a 
bow shock. This shock propagates also in the direction 
perpendicular to the blob's velocity. As the shocked blob 
material heats up it begins to expand backwards and a 
backwards going shock develops as well. The expected 
morphology of the remnant will therefore be of two 
expanding shells which will eventually join, producing 
yet another shock. At late times the shells merge and 
become a single spherical shell. 

Eig. 2 depicts the expected schematic hydrodynamic 
behaviour after the Newtonain transition. This is indeed 
confirmed in the computation. Eig. 4, depicts the density 
contours along the evolution. We observe the expected 
evolution: from two individual blobs via a peanut shape 
configuration with a shock along the equator towards a 
more and more spherical configuration at late times. 



The ratio Zmax/rxy can be approximated by a power 
law as shown in Eig 5. In our simulation this ratio is 
always between 1 and 2 so that the power law fit is very 
inaccurate. This ratio decreases in time as a power law 
with an exponent of —0. 15. Extrapolating this power law 
we see that this ratio reaches a value of 1 at ~ 2. 1 x 10^. 
At this time the shock has a spherical shape with z = 
r^y ^ 15(£5i/^)^/^pc. Even then the shock will not be 
completely spherically symmetric as there would still be 
a ring of shocks around the "equator" where the shells 
have collided. 

Eigures 6 and 7 depict the images of the remnant as 
a function of time and angles of inclination. We show 
images due to bremsstrahlung emission and synchrotron 
emission. The images are constructed assuming that all 
the gas is optically thin in the relevant frequencies. The 
bremsstrahlung luminosity (Eig. 6) was calculated as- 
suming that the volume emissivity is proportional to 
p2gi/2 calculating the synchrotron emissivity 

(Eig. 7) we assumed that both the magnetic field energy 
density and the number density of the relativistic elec- 
trons are proportional to the internal energy density of 
the gas with constant proportionality factors 85 and £e 




FIGURE 5. The ratio between the radius rxy of the shock and 
the z position of the shock. The soUd Une is the best fit power 
law . The ratio will reach a value of 1 at ^ ^ 2. 1 x 10^ yr 



respectively. We further assume that the relativistic elec- 
tron number density is a power law in energy. Under 
these assumptions the volume emissivity is proportional 
to p^8^ [e.g. 20]. In the late images there are two bright 
circles at the lines where the colliding blobs form a hot 
shocked region. In figures 8 and 9 we show the charac- 
teristic emission frequencies. For bremsstrahlung this is 
kT/h where T is the temperature of the gas. For syn- 
chrotron emission we assume 85 =0.1. The characteris- 
tic frequency in this case is the Larmor frequency eB/rUe 
where e, B and rUe are the electron charge, the magnetic 
field and the electron mass respectively. 
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FIGURE 6. Images of the remnant, bremsstrahlung emis- 
sion. The number above each image is the angle of inclination 
in degrees. The images are shown at the same fd as the last 8 
panels of figure 4. 



DISCUSSION 

The long time shape of a GRB remnant is insensitive 
to the exact initial morphology, angular width and den- 
sity of the ejecta. Initially the remnant is highly non- 
spherical. It becomes spherical as time advances and 
the ratio between its height and radius approaches unity 
when 

A/^2.1xl0^ (2) 

This corresponds to 

t^\(fyxEl{\-"\ (3) 

and 

R^n^c{E,,/n)'IK (4) 

After this time it will be difficult to distinguish a GRB 
remnant from a SNR on the basis of its morphology 
alone. 

Using as the observed GRB rate Rgrb = 
10~^yr~^ gal~^ [3] the expected number of non- 



spherical GRBRs per galaxy is: 

0.5yr-(A)-i(^)£V3„-i/3. (5) 

This value depends of course critically on the typical 
beaming factor, It should be compared with the ex- 
pectation of 100 similar aged (10^ yrs) SNRs per galaxy. 
We would expect 20 non spherical GRBRs up to a dis- 
tance of 10 Mpc. The angular sizes of these GRBRs 
would be around a ^arcsec. 



Implications to DEM L 316 

DEM L316 [21] in the LMC looks like two col- 
liding bubbles (see Fig. 10) . It is thought to result 
from a collision between two SNRs. This requires, of 
course, an unlikely coincidence in the timing and the lo- 
cation of the two SNes. An interesting possibility is that 
DEM L 316 is a GRBR. Does this fit our model? DEM 
L316 is far from spherical and has a distinct double 
shell morphology, most similar to our results at ~ 10^ 
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FIGURE 7. Images of the remnant, synchrotron emission. 
The fi are the same as the last 8 panels in figure 6 




FIGURE 8. The characteristic synchrotron frequency as a 
function of fd. 



(see Fig. 4). The /i ratio measured for DEM L316 is 
^1 X 10^. However, according to our results this is far 
after the spherical transition. A GRB remnant would al- 
ready be spherical at this stage. This discrepancy rules 
out the possibility of fitting DEM L 316 with our model 
for a GRB remnant. 
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FIGURE 9. The characteristic bremsstrahlung frequency as 
a function of fd. 




FIGURE 10. DEM L 316 in the LMC (from [21]). 



An Underlying Supernova 

Our model should be modified if the GRB relativistic 
beams are accompanied by an underlying spherical su- 
pernova, as would be expected in the Collapsar model 
[14, 15, 16]. In this case a spherical shell of ~ lOm©, 
the supernova ejecta, will accompany the GRB beams. 
This ejecta propagates at a much lower, Newtonain ve- 
locity, with initial values of ~ lO^km/sec. However, it 
will not slow down while the outer GRB ejecta is piling 
up the external matter and is slowed down. Eventually 
it will catch the GRB ejecta. It is clear that at this stage 
the SN shell would tend to make the GRBR bow shock 
more spherical. To determine when this will happen one 
needs another set of numerical simulations. These are in 
progress now. However, we can attempt to estimate when 
the slower SN ejecta will catch up the slowing down 



GRBR remnant. Assuming that the SN ejecta does not 
slow down (as the GRBR ejecta clears the surrounding 
ISM matter) we find that this will happen at: ^ ^ 3000, 
namely at t ^ 150yr(£5i/^)i/^and R ^ 4pc(£5i/^)i/l 
This happens around the time that the two shells collide 
on the equator. We expect to see an enhanced emission 
due to this collision and then the system will become 
quickly spherical. The expected number of non spherical 
GRBRs and their corresponding sizes would be smaller 
by a factor of 10 then the values estimated earlier for the 
simple evolution of the beamed GRB ejecta. Thus we ex- 
pect one or two non spherical GRBRs with distances up 
to 20 Mpc and their sizes would be around 0.1 ^arcsec. 
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